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We study the parity effect of the particle number in the interference fringes of a Bose-Einstein 
condensate released from a double- well potential. For a coherently splitting condensate in the double- 
well potential, with a decoupled two-mode Bose-Hubbard model, there is well-known phase diffusion 
because of interatomic interactions. After a specific holding time of the double-well potential, the 
phase diffusion will make the interference patterns in the density distribution depend strongly on 
the parity of the total particle number by further overlapping two condensates. This parity effect 
originates from the quantized relative phase about the total particle number. The experimental 
scheme to observe this "even-odd "effect of the particle number is discussed. 



I. INTRODUCTION 

For most classical systems with a large number of par- 
ticles, adding or decreasing one particle will not change 
the fundamental properties of the system. For the quan- 
tum system of dilute Bose or Fermi gases in a harmonic 
trap with a large number of particles, it is also widely 
believed that adding or decreasing one particle will not 
change the thermodynamic and dynamic properties of 
the system. For example, for a Bose-Einstein condensate 
(BEC) described by the Gross-Pitaevskii equation within 
the mean-field model, adding or decreasing one particle 
plays a negligible role in the coupling parameter (which 
is proportional to the overall particle number) for the or- 
der parameter. Recently, it is found theoretically that for 
specific cold atomic system beyond the mean-field model, 
there is a parity effect of the particle number which means 
that adding or decreasing one particle would has different 
quantum effect. The parity effect has been studied theo- 
retically for the quantum decay of Josephson n states [l[ , 
the Berry phase of two-species Bose-Einstein condensates 
(BECs) [2], the tunnel splitting in the Josephson model 
Q, the macroscopic superpositions of phase states 
and more recently in the mesoscopic quantum switching 

To our best knowledge, the parity effect of particle 
number has not been observed experimentally for cold 
atoms, although the parity-dependent tunneling splitting 
was observed for magnetic molecular clusters [y, 0, M, 
9]. ft is still an open question to observe the parity- 
dependent effect for cold atomic system. For example, it 
is very difficult for an experiment to observe the parity- 
dependent effect in the tunneling splitting, because the 
tunneling splitting is extremely small [5|. 

The quantum interference between two condensates in 
a double-well potential pH [Q, O Q OS M, G3, M, 



EH HH HH, H2| has renewed interest in the experimental 
and theoretical studies of the macroscopic quantum co- 
herence effect. In the last few years, there arc significant 
advances in the nonlinear self-trapping of weakly coupled 
condensates [IH [24| , long phase coherence time for two 
separated condensates [25| etc. Stimulated by the re- 
markable experimental advances of BECs in double-well 
potential, we find strong parity-dependent interference 
fringes for ultracold bosonic gases released from a double- 
well potential, based on the phase diffusion in two-mode 
Bose-Hubbard model. For specific holding time of the 
double-well potential to create the phase diffusion, we 
find that the center in the density distribution is a dip 
for even particle number, while it is a peak for odd par- 
ticle number. This provides a way to display directly the 
quantized characteristic of the particle number. 

The paper is organized as follows. In Sec. II, we con- 
sider the phase diffusion of a condensate in a double- 
well potential, with the decoupled Bose-Hubbard Hamil- 
tonian. The parity effect of the particle number is dis- 
cussed based on the coherence property. In Sec. Ill, 
we predict strong parity effect of the particle number in 
the interference patterns of the density distribution. In 
Sec. IV, the influence of the parity effect due to asymme- 
try fluctuations of the double-well potential and particle 
number squeezing is studied. In the last section, we give 
a brief summary and discuss the application of Feshbach 
resonance to observe more clearly the parity effect. 



II. PHASE DIFFUSION AND COHERENCE 
PROPERTY OF A CONDENSATE IN A 
DOUBLE- WELL POTENTIAL 

For a condensate in a double-well potential with neg- 
ligible tunneling, the system can be described by the fol- 
lowing decoupled Bose-Hubbard Hamiltonian. 
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Here = a\a,i is the particle number operator for the ith 
condensate, a, and a\ are bosonic annihilation and cre- 
ation operators for the zth condensate. £1 — 62 represents 
the bias potential between two sites, while U denotes col- 
lisional interaction energy. The initial quantum state in 
this double-well potential is assumed as 
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The evolution of the quantum state is then 
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We see that there is a phase diffusion [2J, |27|, |2J, |29| 
in the quantum state (th)) because of the exponen- 
tial factor. In addition, from the exponential factor in 
I* (</,)), the time evolution of |\fr (t^)) displays a peri- 
odic behavior of T = Anh/U . To give a clear discussion 
of the parity effect, in getting Eq. (|2|), we consider the 
case of e\ = e-2- 

The above evolution of |\& (th)) may be realized by 
first preparing two coherently separated condensates in a 
double-well potential with large Josephson tunneling, so 
that the quantum state can be described well by Eq. (2]| . 
The central barrier is then increased non-adiabatically so 
that the tunneling between two condensates can be omit- 
ted, while the quantum state of the system still takes the 
form given by Eq. Ip]). This is similar to the experimen- 
tal studies of the collapse and revival of the quantum 
coherence of the matter wave packet in an optical lattice 
[30l | . This scheme was also applied in Ref . Q to consider 
the macroscopic superpositions of phase states. 

A direct way to display the coherence property of 
I "J" (th)) is to consider the time evolution of K (th) — 

(* (th)\ fai«2 + a^cixj |\P (th)) /N, which reflects the 

phase coherence and interference effect between two con- 
densates. From the expression ([3]), we have 
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The above equation can be rewritten as 
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Based on the binomial theorem, we have 
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FIG. 1: Shown is the time evolution of the phase coherence 
factor K (tho)- Here t^o is in unit of T in the text. Besides 
the periodic behavior of K (tm), we see significantly different 
behavior in K (tho) for odd and even particle number. 



In this situation, we have the following simple expression 

(cos (Uth/h))"- 1 
Ku(th) = o ■ ( 7 ) 



From the above results, we have 



K (t h0 ) = (cos (intho)) 



N-l 
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where tho = th/T. For large N, by using the identity 
limit lirri/v^oo (1 — x/N) — e~ x , the above expression 
can be approximated very well as 

n>0 
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where n is an integer. In this situation, the width of the 
peaks or dips in K (tho) is about 1/2-k^/N. At tho = 0, 
there is ideal phase coherence between two condensates. 
As shown in Fig. 1(a) for N = 101, with time in- 
creasing, the phase coherence disappears very rapidly for 
tho > l/2w\/N. With further time increasing, however, 
there is a revival of the phase coherence. This behavior 
is understood from the periodic behavior of the quantum 
state \$(t h0 )). 

A unique parity effect of the particle number is shown 
by the expression (8| of K (tho)- For specific holding time 
such that cos (Ant ho) = — 1, we see that K (tho) = 1 for 
odd particle number, while K (tho) = — 1 for even particle 
number. This parity effect is shown further in Fig. 1(a) 
and Fig. 1(b) for N = 101 and N — 100, respectively. 



III. PARITY EFFECT OF THE PARTICLE 
NUMBER IN THE DENSITY DISTRIBUTION 



The above parity effect can be considered further 
by overlapping two condensates after switching off the 
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FIG. 2: (Color online) Shown is the density distribution for 
different holding time t^o and fixed time of flight tfo = 1. The 
solid (dotted) line is the density distribution of the particle 
number N = 100 (N = 101). Clear parity effect of the particle 
number in the interference patterns of the density distribution 
is shown in Figs. 2(c)-2(e). 



double- well potential. After a holding time of tuo, the 
quantum state is I'F (tho))- Then, we calculate the den- 
sity distribution after a time of flight tf by switch- 
ing off the double- well potential. In our calculations, 
the initial wave functions of two condensates are as- 
sumed as i Pl (as) = e -(^+d/2) 2 /2 CT 2 / 7r i/4 (J i/2 and V2 ( x ) = 

e -(^- d /2) /2a ^1/4^1/2^ respectively. In this situation, 
the distance between two condensates is d. In the follow- 
ing calculations, we will adopt the length unit d, energy 
unit Ed = h 2 /2md 2 and time unit Td = h/Ed- After a 
time of flight tf = i/o^d, the density distribution of the 
system is given by 

N r 2 2 

n(x,tho,tfo) = — \\fa (x,tf )\ + \fa(x,t f0 )\ + 
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For N = 100 and N = 101, in Figs. 2(a)-(f), we give the 
density distribution for fixed time of flight tto = 1 and 
different holding time tho- In our calculations, a = 0.1. 
The factor (cos (intho)) N 1 in the interference term of 
the above equation gives several unique results: 

(i) There are collapse and revival of the interference 
fringes in the density distribution. 

(ii) At tho = j/2 — 0.25 (j is a natural number), the 
density distribution depends strongly on the parity of 
the total particle number. As shown in Fig. 2(d), the 
center of the density distribution (solid line) is a dip 
for even particle number, while the center of the den- 
sity distribution (dotted line) is a peak for odd parti- 
cle number. This can be understood further by noticing 
that n = N\<t) 1 - fa\ /2 for even particle number, while 
n = N |0i + fa\ /2 for odd particle number. For both 
situations, the density distribution is symmetric, which 



means a clear parity effect about the total particle num- 
ber m. 

At tho — j/2 — 0.25, the factor [cos (47r^o)] Ar_1 m the 
interference term can be rewritten as e 171 ^ -1 ). In this 
situation, the phase diffusion of the left and right conden- 
sates is zero. However, there is a relative phase between 
two condensates. This relative phase is quantized about 
the total particle number. For even particle number, the 
relative phase is tt, while the relative phase is 2-7T for odd 
particle number. 

Before the atomic cloud is imaged on a CCD camera, 
we do not know the parity of the particle number. At 
tho = i/2 — 0.25, the parity effect would lead to a ran- 
dom behavior in the density distribution. Averaging the 
density distribution in different experiments, there would 
be no interference patterns. However, in a single experi- 
ment, there will be clear interference patterns, and there 
is 50% probability to observe a peak (dip) in the center 
of the density distribution. This is significantly differ- 
ent from the holding time of tho = j/2 — 0.5, where one 
would always observe the same interference patterns in 
the density distribution. 

These analyses at tho = j/2 — 0.25 can be considered 
further by the following density matrix 



p = P(N)\N) (N\ 



(11) 



Here P (N) is the probability to have N atoms in the 
condensate. It can be assumed as a Gaussian distribu- 
tion about the average particle number TV. The above 
expression can be rewritten as 

p= P(N)\N){N\+ ]T P(N)\N){N\. (12) 

N—even N—odd 

The ensemble average of the density distribution is then 



;[x,t ho ,tfo) = Tr 



N r 
~2 



\fa(x,tf )\ +\fa(x,tf )\' 



(13) 



It is clear that there is no interference term in the ensem- 
ble average of the density distribution. In a single-shot 
density distribution, however, there would be interference 
fringes with random relative phase 7r or 2-7r. 



IV. INFLUENCE OF THE PARITY EFFECT 
DUE TO ASYMMETRY OF THE DOUBLE- WELL 
POTENTIAL AND PARTICLE NUMBER 
SQUEEZING 

As shown in preceding section, the density distribu- 
tion depends strongly on the parity of the overall par- 
ticle number N. It is a natural problem whether slight 
asymmetry fluctuations of the double-well potential will 
destroy the parity effect. Generally speaking, the asym- 
metry fluctuations of the double-well potential will lead 
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to a random relative phase and fluctuations of the aver- 
age particle number in two sites. To consider both effects, 
we study the following more general initial state 

|*fo = 0)> = ^2=(aa\ + ^f\Q). (14) 



At tho = j/2 — 0.25, we have 



A' 



K12 = j^J2VW^TTY) 



-[{l-N/2) 2 + (l-N/2-l) 2 ]s 2 /N e i7t{N+l) _ (19) 



The normalization condition requests that \a\ +\/3\ — 2. At this holding time, we can also get 
At a holding time tho, the density distribution after a 
time of flight t fo is given by 



K 2 i = 



n{x,t h0 ,t f0 ) = — ||a| 2 \4> x (x,tf )\ 2 + \(3\ 2 \(f> 2 (x,tf )\' 
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At the holding time tho = j/2 — 0.25, the density dis- 
tribution is then 

n(x,t h0 ,t f0 ) = — ||a| 2 |0i {x,t f o)\ 2 + \(3\ 2 \(f> 2 (x,t f0 )\ 2 

+2 (-if' 1 Re {a* 13^*1 (M/o) <f>2 (*,t/o)]} ■ (16) 

We see that there is still strong parity effect about the 
total particle number when the asymmetry of the av- 
erage particle numbers in two sites is considered. For 
tho = j/2 — 0.25, n = N \a<f>i — fi<\>2\ /2 for even parti- 
cle number, while n = N \auf>i + Pfal /2 for odd particle 
number. 

The factor a* (3 includes a random relative phase if the 
fluctuations of the bias potential E\ — £2 are considered. 
The bias potential £1 — £2 between two sites will gives a 
factor e ~ 4 ( £l ~ £2 ) th /'\ Experimentally, to observe the par- 
ity effect, this requests that the fluctuations of the bias 
potential should satisfy \Ae\th/H « ir. This request 
can be satisfied in the present experimental technique 
HI, HH ■ These analyses show that slight asymmetry 
fluctuations of the double-well potential would not lead 
to serious problem in observing the parity effect. 

In the non-adiabatical increasing of the central bar- 
rier so that the tunneling between two condensates can 
be omitted, it is possible that there would be a parti- 
cle number squeezing. In this situation, we consider the 
following quantum state 



I* (th = 0)> 



2=0 



7e -a-JV/2)^VJV | Z ,jv- . (17) 



Here s is the squeezing factor, and 7 is a normalization 
factor. For this quantum state, we have 
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We see that at tho — j/2 — 0.25, even particle number 
corresponds to a relative phase of w, while odd particle 
number corresponds to a relative phase of 27r. This means 
that the squeezing factor will not destroy the parity effect 
in the interference fringes of the density distribution. Of 
course, there is a request that the quantum state should 
not be a Fork state in the extreme squeezing situation. 



V. SUMMARY AND DISCUSSION 

In summary, a strong parity effect of the total parti- 
cle number in the density distribution of a condensate 
released from a double-well potential is predicted in this 
work. The asymmetry fluctuations of the double-well 
potential and particle number squeezing are discussed, 
and our studies show that they would not lead to serious 
problem in observing the parity effect. The parity effect 
is due to the quantized relative phase about the total par- 
ticle number for specific holding time of the double-well 
potential. Another example of the strong parity effect 
about the total particle number is the exchange effect for 
a complex of JV spin 1/2 fermions [33[ . The parity effect 
of the complex lies in that it is a boson for even JV, while 
it is a fermion for odd JV. 

In this work, similarly to the phase diffusion consid- 
ered in Ref. [27|, we ignore the phase diffusion in the 
increasing of the central barrier by assuming that the 
central barrier is increased fast enough, and we also ig- 
nore the phase diffusion in the time of flight. During 
the time of flight, the interaction energy per particle will 
decrease rapidly, and thus the phase diffusion may be 
ignored. In Ref. the phase diffusion in the time of 
flight is also ignored. Nevertheless, overcoming these two 
phase diffusions would contribute to more clear observa- 
tion of the parity effect. One way to overcome these two 
phase diffusions would be the application of Feshbach 
resonance which has been used in atom interferometry 
[341 [H[ . The scheme with the application of Feshbach 
resonance would be: (i) One first prepares two coherently 
separated condensates in a double-well potential, (ii) Af- 
ter turning the s-wave scattering length almost to zero via 
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a magnetic-field Feshbach resonance, the central barrier 
is increased adiabatically so that the tunneling between 
two sites can be omitted, (iii) Holding the double-well 
potential, one adiabatically increases and then decreases 
the s-wave scattering length almost to zero so that the 
condition f* h ^-jp-dt — n is satisfied, (iv) Switching off 
the double-well potential, the ideal gas is then imaged 
after a time of flight to display the parity effect. 
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